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Abstract. - Competing ferro- and antiferromagnetic exchange interactions may lead to the for- 
mation of bound magnon pairs in the high-field phase of a frustrated quantum magnet. With 
decreasing field, magnon pairs undergo a Bose-condensation prior to the onset of a conventional 
one-magnon instability. We develop an analytical approach to study the zero-temperature prop- 
erties of the magnon-pair condensate, which is a bosonic analog of the BCS superconductors. 
Representation of the condensate wave-function in terms of the coherent bosonic states reveals 
the spin-nematic symmetry of the ground-state and allows one to calculate various static prop- 
erties. Sharp quasiparticle excitations are found in the nematic state with a small finite gap. 
We also predict the existence of a long-range ordered spin-nematic phase in the frustrated chain 
material LiCuV04 at high fields. 



Introduction. — Frustrated spin systems are interest- 
ing in general, as their zero- and low-temperature proper- 
ties are governed by quantum fluctuations. These strong 
fluctuations may inhibit the formation of long-range mag- 
netic order, stabilizing instead a disordered spin liquid. 
Following Anderson [T], much of the interest in the past 
' decades has been focused on the investigation of various 
\ quantum spin-liquid states . Another interesting pos- 
• sibility is the appearance of unconventional magnetic or- 
der characterized by partial breaking of the spin-rotational 
symmetry 0(3). A specific example is provided by spin- 
, nematic states, which are analogous to the ordered phases 
of needle-like molecules in liquid crystals. Spin-nematic 
phases have been discussed phenomenologically in [HIS], 
whereas identification of the relevant microscopic mecha- 
nism remains a challenging theoretical problem. 

It was suggested long time ago [6^ that a sizable bi- 
quadratic exchange (Si • S^)^ in magnetic insulators with 
S* > 1 may stabilize a quadrupolar phase with vanishing 
sublattice magnetization (S) = 0, but a nonzero second- 
rank tensor, e.g., (S*^) 7^ {Sy) = {Si). In most real com- 
pounds the biquadratic exchange is, however, rather small. 
Recently, the interest in the biquadratic mechanism for the 
spin-nematic order has been revived in connection with the 
experiments on cold atomic gases [2l|8! and on the disor- 



dered magnetic material NiGa2S4 [91111). 

In this Letter we explore an alternative mechanism for 
the spin-nematic ordering based on competition between 
ferro- and antiferromagnetic interactions in magnetic in- 
sulators with an arbitrary value of the local spin includ- 
ing S — 1/2. Specifically, the mechanism operates in 
strong magnetic field and is based on the formation of 
bound magnon pairs in the fully polarized state |12fll5| . 
This scenario has been studied numerically in a num- 
ber of works on the so-called ferromagnetic J1-J2 chain 
model p^H20] and its generalization to two dimensions 
(2D) [2TJ[22]. Clear numerical evidence was found for the 
critical quadrupolar correlations below the saturation field 
for the ID case. Note, however, that there is no true long- 
range order, nematic and otherwise, in one-dimensional 
quantum magnets at zero temperature. Therefore, a num- 
ber of important questions on stability of the ordered ne- 
matic state and its excitation spectra are not answered by 
studying the purely ID model. 

The purpose of this Letter is to develop a simple ana- 
lytical framework to treat the ground-state properties and 
low-energy magnetic excitations in the phase with a long- 
range spin-nematic order. Our description of the conden- 
sate of bound magnon pairs resembles in many aspects 
the BCS theory for the condensate of bound electron pairs 
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Fig. 1: (Color online) Energy-field diagram for a frustrated 
quantum magnet close to the saturation field. Dot-dashed lines 
show lowest one- and two-magnon states. Solid lines represent 
the ground state energy for the one-magnon (spin-cone) and 
the two-magnon (spin-nematic) condensate. 

in superconductors. In addition, we predict that a spin- 
nematic phase must exist at high fields in the frustrated 
chain material LiCuV04 p5H^ . 

The phenomenon of magnon pair condensation has a 
close relationship to the old problem of particle versus 
pair-superfluidity in an attractive Bose gas [52]. A com- 
mon outcome is the density collapse prior to the pair- 
condensation transition [30]. Unrestricted growth of the 
local magnon density in spin- 1/2 antiferromagnets is cured 
by their hard-core repulsion. In addition, reduced dimen- 
sionality of a spin subsystem found in many real materials 
helps to stabilize bound pairs and creates favorable con- 
ditions for their condensation. Experiments on LiCuV04 
and other related compounds in high magnetic fields may, 
therefore, lead to the first observation of such an exotic 
off-diagonal long-range order in solid-state systems. 

In order to demonstrate the occurrence of magnon pair 
condensation at high magnetic fields we start with a gen- 
eral quantum Heisenberg antiferromagnet on an A'^-site 
lattice 

i,r i 

where r = — r^. In the following S ~ 1/2 is set for 
definiteness. In strong magnetic fields the Zeeman energy 
dominates over the exchange interactions and stabilizes 
the fully polarized state |0) = |ttt- ■ •)• This state is the 
vacuum for single spin-flips or magnons 

|ii) = ;^E^"'''""'^rlo) (2) 

with the excitation energy 

r 

where Jq = X^r "^('")^*'*'"- -'-'^ ordinary antiferromagnets 
spin-flips repel each other. Then, once the band gap in 



vanishes at a certain Q, Jq = min{Jq}, an antiferromag- 
net undergoes a second-order transition into a canted spin 
structure at the saturation field 

^.1 = ^ (-^0 - Jq) , (4) 

as illustrated in Fig. [T] The antiferromagnetic state be- 
low Hsi can be regarded as a Bose-condensate of single 
magnons [5T[|^ . 

Two-magnon bound states. — The conventional 
scenario for an antiferromagnetic transition in a strong 
magnetic field may change if some of the exchange bonds 
are ferromagnetic. In this case two spin-flips occupying the 
same bond with J(r) < lower their interaction energy 
and may form a bound pair [12l[T3]. To treat the bound 
state problem we follow the standard approach [331135] and 
define a general two-magnon state 

with /,y — fji being the magnon pair wave-function. Sep- 
arating the center of mass motion fij ~ e'''('''+'"j^/^ /k(r) 
and calculating the matrix elements of H for states (O we 
obtain the Bethe-Salpeter equation 

(£2 - ek/2+q- ek/2-q) /k(q) = (6) 

^ 2/V E("^P+q+ "^P-q^ Jk/2+q- •4/2-q)/k(p), 
P 

where /k(q) is the Fourier transform of /k(r) and £2 mea- 
sures the magnon pair energy relative to the energy of 
the ferromagnetically polarized state. The above equa- 
tion extends the previous theories [12l[l3l[T5] to an arbi- 
trary geometry of exchange interactions and with a trivial 
replacement of Eq it also remains valid for an arbitrary 
value of spin S. 

While the subsequent theoretical consideration is en- 
tirely general, we introduce now for illustration a specific 
spin model shown in Fig. [21 which is related to the quasi- 
ID helical antiferromagnet LiCuV04. This material con- 
sists of planar arrays of spin-1/2 copper chains with a fer- 
romagnetic nearest-neighbor exchange Ji < and an an- 
tiferromagnetic second- neighbor coupling J2 > 0. Chains 
are linked by diagonal bonds, whereas interplanar interac- 
tions are an order of magnitude smaller. Neutron scat- 
tering measurements provide the following estimate for 
the exchange parameters in LiCuV04: Ji = —1.6 nieV, 
J2 = 3.8 meV, and J3 = —0.4 meV [24]. The importance 
of quantum effects in this material is revealed by a small 
value of ordered moments ~ 0.3/^s in zero magnetic field 
[23] . The purely ID J1-J2 model has been studied inten- 
sively in the past [T2 l[T6H20] though no results exist for a 
realistic planar model. 

To solve the integral equation ([6]) we expand /k(q) into 
lattice harmonics and obtain a finite algebraic system. 
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Fig. 2: (Color online) Two dimensional array of copper ions in 
LiCuV04 with principal exchange couplings. 



Bound states exist for O.Stt < ky < n with the mini- 
mum of e2(k) at K = (7i',7r). The binding energy de- 
fined by e2(K) = 2eQ — Eb is found numerically to be 
Eb ~ O.O3OJ2. In the absence of bound states the two- 
magnon continuum has a gap that is twice larger than the 
lowest one-magnon energy eq = H — Hgi- The two gaps 
vanish, therefore, in the same magnetic field. When bound 
states are present, condensation of magnon pairs starts in 
a higher magnetic field: 



Hsi + - Eb 



(7) 



see Fig. [TJ The condensation field for LiCuV04 is cal- 
culated to be Hs2 — 47.1 T (5 = 2), whereas the single 
magnon branch softens at Hsi = 46.5 T. The relation 
Hs2 > Hsi holds up to J3 « —0.6 meV. Hence, the con- 
clusion about the magnon pair condensation in LiCuV04 
is rather robust and should not be affected by a possible 
uncertainty in the experimental coupling constants ^23]. 

We finish our analysis of the linear problem by present- 
ing the wave-function of the lowest energy bound pairs in 
the momentum representation: 



/K(q) 



XJ2 cos ( 



2^2(1 — cos2qy) + 4J3 sinqx sinqy + Se ' 



(8) 



where a numerical constant 6e — O.I3OJ2 is related to the 
binding energy and A is the normalization factor. In real 
space the function e*^'"/^/K(i") has the odd parity under 
the reflection y — > —y and vanishes, therefore, at r = 0. 

Coherent condensate of magnon pairs. Below 
Hs2 the bound magnon pairs acquire negative energy and 
start to condense. It is convenient at this point to trans- 
form from spin- 1/2 operators to the Holstein-Primakoff 
bosons: 



5f = - - aja. 



(9) 



expanding subsequently square-roots to the first order in 



The many-body state with a macroscopic number of 
the lowest-energy pairs below Hs2 can be expressed as the 
coherent boson state of the pair creation operator: 



|A) 



-A'|A|^/2 



exp 



aja] 



|0) . (10) 



1,3 



Here fij — 6*^'^'"'+'"^^/^ /K(r) is the wave-function of the 
lowest energy pairs and A is the complex amplitude of 
the condensate. The state ([T0| is a bosonic equivalent 
of the BCS pairing wave-function for fermions and must 
be regarded as a variational ansatz, which can be further 
improved by taking into account pair-pair correlations. 

We use the ground-state wave-function pUj) to compute 
simple boson averages: 

/ \ _ n / \ _ A/K(q) 

[aq) - U , ^aK/2+qaK/2-q; - Y^^~\A^J^{^ ' 

nK/2+q = (aK/2+qaK/2+q) = T"^]Xp^(q) '"^^^^ 

as well as a more complicated four-boson correlator: 

("p/2 + q«p/2-q«P/2+q'«p/2-q') = |AP<5p,K 

^ n _L ^ ^ /K(q)/K(qO 

X (1 + nK/2+q + nK/2+qO i_|A|4;^(q);^(q,) 



+ 1^1" iKq' + (1 + "p/2+q + Wp/2-q) 



2 f p-K 



2 ( P-Vi 



l-|A|V^(^+q)/i(^-q) 



(12) 



From these one can derive various spin correlators. In 
particular, the absence of the single-magnon condensate 



(flq) = translates into 







(13) 



which signifies a lack of a usual antifcrromagnetic order 
parameter. The transverse and longitudinal spin correla- 
tions are given in the leading order by 



lAI 



I 



f^i,, {SS!SS])^\A\'\n,f (14) 



with SS^ = — {Sf). In accordance with the behavior 
of the bound-state wave- function /K(r) the two correla- 
tors decay exponentially as | r^ — r j | — ^ cxo . The transverse 
magnetic structure factor 5^ (q) has no Bragg peaks and 
exhibits diffuse liquid-like spin correlations with a char- 
acteristic shape in momentum space determined by the 
bound-state wave-function: 



lAI 



/K(f + q) + /l:(f-q) 



(15) 



The longitudinal response S'^^(q) being formally of the 
same order in |A| appears to be much weaker than >5'^(q) 
in the vicinity of the saturation field Hg2 ■ 
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The general definition of the spin nematic order param- 
eter in the 0(2)-symnietric case is 

Qtf = \ {S-S^ + S^Sf) ~ \ 5^1, {Si- ■ Sf) (16) 

where i,j belong to a nearest-neighbor bond and a, (3 — 
X, y. The quadrupolar tensor Qf^^ acquires a nonzero ex- 
pectation value in the presence of the pair condensate: 

Q^r+'Q^'-l(S+Sp^^f,,. (17) 

The phase of the condensate amplitude A determines the 
orientation of the spin-nematic director in the x-y plane. 
The director forms a periodic structure in the real space 
determined by the momentum K. 

Using explicit expressions for the spin correlators, one 
can calculate the ground-state energy. Below, we shall 
use for this purpose a complementary approach, which 
is analogous to the Bogoliubov method in the theory of 
superconductivity and yields the spectrum of quasiparticle 
excitations simultaneously with the static properties. 

Mean-field approach. — The bosonic equivalent of 
the spin Hamiltonian is obtained upon substitution of Q 
into (dJ. We restrict ourselves to quadratic and quartic 
terms in a^'s and define two types of mean- field averages 
for each exchange bond: 

Ay = {aiaj) , nr = (aja^) , (18) 

and the magnon density n = {a\ai). The anomalous cor- 
relator is further factorized as A^j = e*^'''*+'"j'/^Ar. Both 
Ar and rir are even real functions of r with a proper 
choice of gauge. Performing the mean-field decoupling in 
the interaction term we obtain a quadratic form, which is 
then diagonalized with the canonical transformation. This 
yields the energy of one-magnon excitations 

£K/2+q = fci-^ J - n - n,,) sin ^Kr sinqr, 

r 

eq-/42~~B2, Bq-^ J(r)ArCOsqr , (19) 

r 

Ac^ = H — ^ J{^)i\ — n — nr)(l — cos ^Kr cosqr). 

r 

In accordance with the exponential decay of spin correla- 
tions (jl4p . the excitation spectrum acquires a gap in the 
presence of the magnon pair condensate. The above ex- 
pressions are used to calculate bosonic averages and to 
obtain a closed form of the self-consistent equations: 

^"■ = E^^o^q^' "r = ^^cos(iK-)-q)r. (20) 

Z6q Z6q 

In the limit H — > Hs2 one finds Ar ^ rir ^ A^, while 
the linearized equation for Ar transforms directly into the 
bound-state equation ([6]). 
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Fig. 3: (Color online) Dispersion of a single-magnon branch 
in LiCuV04 in the fully polarized state {Hs2 ~ 47 T) and in 
the state with the magnon pair condensate. Field values for 
different curves from top to bottom are listed on the plot. The 
inset shows the vicinity of the magnon gap. 

We have solved self-consistently the set of equations 
and (pOj) and calculated the ground-state energy for 
the spin model of Fig. [2] assuming the same symmetry of 
bond variables Ar and as in the magnon pair coherent 
state pil)l . With decreasing external field, the pairs over- 
lap more appreciably and at a certain point give way to 
a conventional one-particle condensation. Comparing the 
ground-state energy of the pair condensate with the en- 
ergy of a simple spin-cone structure we find the first-order 
transition at He ~ 44.5 T as illustrated schematically in 
Fig. [T] The spin-nematic state has the lowest energy in a 
finite range of fields He < H < Hs2, which extends well 
below the condensation field Hgi for single magnons. 

The ground-state energy calculation allows to determine 
the slope of the magnetization curve M (H) ai He < H < 
Hg2- In ordinary quasi-lD antifcrromagnets, M{H) devi- 
ates from a straight line as H ^ Hg resembling the square- 
root singularity of a single quantum spin chain. Our 
mean-field calculations for the high-field nematic phase in 
LiCuV04 yield instead the slope dM/dH « 0.38Msat/J2, 
which amounts to only 54% of the slope of the classical 
magnetization curve. The quantum corrections beyond 
the mean-field approximation should somewhat modify 
this value. However, we expect them to be small for the 
same reason as the suppression of critical fiuctuations in 
the BCS superconductors. Indeed, the size of the bound 
magnon pairs is rather large extending to ^ ~ 10 inter- 
atomic spacing in the direction of chains. Already for 
small magnon densities each bound pair is surrounded by 
many neighboring magnon pairs, which enforces the mean- 
field behavior. Therefore, a distinct signature of the high- 
field nematic phase in LiCuV04 will be a sharp change in 
the slope of the magnetization curve. 

The dispersion of one-magnon excitations found to- 
gether with the ground-state energy in the self-consistent 
calculation is presented in Fig. [3] In the fully polarized 
state at H > Hs2 ~ 47 T, variation of the applied field 
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results in an overall shift of the magnon energy accord- 
ing to Eq. At H = Hs2 magnons have a small gap 
Ag = Hs2 — Hsi fa 0.06 meV. The field dependence of 
Ek changes drastically in the presence of the magnon pair 
condensate. Decreasing the field modifies the shape of Ek 
but the excitation gap remains practically unchanged, 
see the inset in Fig. [3] The lowest field H = 40 T used in 
Fig. [3] is below the transition field He into the spin-cone 
magnetic structure. The corresponding dispersion curve 
Ek illustrates that the spin-nematic state remains locally 
stable even below He- This is in contrast with the previous 
scenario suggested for an attractive Bose gas, for which the 
pair condensate was assumed to become unstable due to 
a softening of the single-particle branch [291130] . 

The motion of the spin-nematic order parameter pro- 
vides an additional gapless branch of collective excitations, 
which should yield a nonzero dynamical signal in the lon- 
gitudinal channel. Transverse spin-spin correlations in the 
nematic phase are dominated by unpaired magnons (jl9p . 

Finally, let us also comment on a low-field state at 
H < He- The considered scenario of a direct transi- 
tion between the spin-nematic state and the conventional 
canted antiferromagnetic phase applies most certainly to 
the 2D model [21], which exhibits a magnetic ordering 
at Q = (7r,0) in zero field. In the case of weakly cou- 
pled chains, the low-field phase might be more complicated 
than the simple conical spin-structure used above as an ex- 
ample. The NMR measurements [26lli28] indicate that the 
intermediate-field phase, which is observed in LiCuV04 
above 8 T [3^, has predominant longitudinal SDW-type 
correlations between local spins. This experimental find- 
ing agrees, in principle, with the numerical results for a 
single spin chain [T7^I19| . Understanding the fate of such 
a ID phase in the presence of interchain couplings as well 
as its relation to the transverse nematic order remains an 
open theoretical problem. 

To summarize, we have presented the analytical descrip- 
tion for the Bose-condensate of bound magnon pairs in a 
frustrated quantum magnet in high magnetic fields. The 
theory applies to a number of real magnetic compounds 
with competing ferro- and antiferromagnetic interactions. 
After the first version of the present work has appeared, 
we learned about the experimental observation of a new 
phase in LiCuV04 in the field range 41-44 T [36]. 
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